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Abstract 

We consider a class of nonlinear Schrodinger equation in two space dimensions with an attractive potential. 
The nonlinearity is local but rather general encompassing for the first time both subcritical and supercrit- 
ical (in L 2 ) nonlinearities. We study the asymptotic stability of the nonlinear bound states, i.e. periodic 
in time localized in space solutions. Our result shows that all solutions with small initial data, converge 
to a nonlinear bound state. Therefore, the nonlinear bound states are asymptotically stable. The proof 
hinges on dispersive estimates that we obtain for the time dependent, Hamiltonian, linearized dynamics 
around a careful chosen one parameter family of bound states that "shadows" the nonlinear evolution of 
the system. Due to the generality of the methods we develop we expect them to extend to the case of 
perturbations of large bound states and to other nonlinear dispersive wave type equations. 

1 Introduction 

In this paper we study the long time behavior of solutions of the nonlinear Schrodinger equation (NLS) with 
potential in two space dimensions (2-d): 

id t u(t,x) = [-A x + V(x)]u + g(u), tel, xeM 2 (1.1) 
u(0,x) = u (x) (1.2) 

where the local nonlinearity is constructed from the real valued, odd, C 2 function g : K i— > M satisfying 

5 (0) = g'(0) = and \g"(s)\ < C(\s\ ai + |sp), a e K, \ < «i < a 2 < oo (1.3) 
which is then extended to a complex function via the gauge symmetry: 

g{e w s) = e w g(s), 9 G K. (1.4) 

The equation has important applications in statistical physics, optics and water waves. It describes certain 
limiting behavior of Bose-Einstein condensates [8j 115] and propagation of time harmonic waves in wave guides 
[HI [Tni HH] • In the latter, t plays the role of the coordinate along the axis of symmetry of the wave guide. 

It is well known that this nonlinear equation admits periodic in time, localized in space solutions (bound 
states or solitary waves). They can be obtained via both variational techniques [TJ (55J HI] an d bifurcation 
methods [501 (55J [T3]- Moreover the set of periodic solutions can be organized as a C 2 manifold (center 
manifold), see [HI [12] or next section. Orbital stability of solitary waves, i.e. stability modulo the group of 
symmetries u i— > e~ te u, was first proved in [551 13H], see also p. [TUl[54] . 

The main result of this paper is that solutions of (|l.ip - (|1.2j) with small initial data asymptotically converge 
to a bound state, see Theorem 13.11 While asymptotic stability results for bound states in NLS have first 
appeared in the work of A. Soffer and M. I. Weinstein [551 HSj, and continued in [5D1 [551 [5[ [3] HI [71 HI], 
our main contribution is to allow for subcritical and critical (L 2 ) nonlinearities, \ < ol\ < 1 in (|1.3p . To 
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accomplish this we carefully project the nonlinear dynamics onto the center manifold of bound states and 
use linearization around this time changing projection to study the motion in the radiative directions, i.e. 
directions that are not in the tangent space of the center manifold. Previously, linearization around a fixed 
bound state has been used, see the papers cited above. By continuously adapting the linear dynamics to 
the actual nonlinear evolution of the solution we can more precisely capture the effective potential induced 
by the nonlinearity g into a time dependent linear operator. Once we have a good understanding of this 
time dependent linear dynamics, i.e. we have good dispersive estimates for its semigroup of operators, see 
Section 2] we obtain information for the nonlinear dynamics via Duhamel formula and contraction principles 
for integral equations, see Section [3] Note that we have recently used a similar technique to show that in the 
critical (cubic) case, with g(s) = s 3 , s <E K, the center manifold of bound states is an attractor for small 
initial data, see [13j . In this paper the technique is much refined, we use a better projection of the dynamics 
on the center manifold and sharper estimates for the linear dynamics. The refinements not only allow us to 
treat a much larger spectrum of nonlinearities including, for the first time, the subcritical ones but also allow 
us to obtain actual convergence of the solution to a bound state. 

However, the main challenge for our approach is to obtain good dispersive estimates for the semigroup of 
operators generated by the time dependent linearization that we use. This is accomplished in Section |4] via 
a perturbative method similar to the one we developed in [13] . As described in that section, we could have 
obtained sharper estimates by using a generalized Fourier multiplier technique to remove the singularity of 

lle^-^IU^ ~ It]" 1 , 

see [12] Section 4]. We chose not to do it because it requires stronger hypotheses on V without allowing us 
to enlarge the spectrum of nonlinearities that we can treat. 

Finally, we remark that our method is quite general, based solely on linearization around nonlinear bound 
states and estimates for integral operators with dispersive kernels. Therefore we expect it to generalize to the 
case of large nonlinear 2D ground states, see for example [7], the presence of multiple families of bound states, 
see for example [27], or to the case of time dependent nonlinearity, see [6]. In all three cases our method 
will not only allow to treat the less dispersive environment, 2D compared to 3D, but it should greatly reduce 
the restrictions on the nonlinearity. The first author and collaborators are currently working on adapting the 
method to other dimensions and other dispersive wave type equations. The work in 3-D is complete, see [12] . 

Notations: H = -A + V; 

L' p = {/ : K 2 i— > C | / measurable and J R2 \f(x)\ p dx < oo}, 1 < p < oo, endowed with the standard norm 

WIWlp = (J* R 2 \ f(x)\ p dx) 1 ^ P ) while for p = oo, L°° = {/ : M 2 i— > C | / measurable and essup|/(x)| < oo}, and 
it is endowed with the norm: = es sup|/(x)|; 

< x >= (1 + l^l 2 ) 1 / 2 , and for a € M, 1 < p < oo, denotes the L p space with weight < x > P<T , i.e. 
the space of functions f(x) such that (< x > a f(x)) p are integrable endowed with the norm ||/(:e)||i,»> = 
|| < x > a f(x)\\ p , while for p = oo, denotes the vector space of measurable functions f(x) such that 
essup| < x > a f{x)\ < oo endowed with the norm ||/(:r)IU?= — II < x >a f{ x )\\L°°', 

(/; <?) — Jr2 f{x)g{x)dx is the scalar product in L 2 where ~z = the complex conjugate of the complex 
number /; 

P c is the projection associated to the continuous spectrum of the self adjoint operator H on L 2 , rangeP c = 

7~Lo] 

H n denote the Sobolev spaces of measurable functions having all distributional partial derivatives up to 
order n in L 2 , \ \ ■ \\un denotes the standard norm in this spaces. 

2 Preliminaries. The center manifold. 

The center manifold is formed by the collection of periodic solutions for (|1.1|) : 

u E (t,x) = e- lEt ^ E {x) (2.1) 
where E £ R and ^ ipE & H 2 (M. 2 ) satisfy the time independent equation: 

[-A + V]^e + 9(1>b) - EipE (2.2) 
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Clearly the function constantly equal to zero is a solution of (|2.2[) but (iii) in the following hypotheses on the 
potential V allows for a bifurcation with a nontrivial, one (complex) parameter family of solutions: 



(HI) Assume that 

(i) There exists C > and p > 3 such that: 

|V(a;)| < C < x >~ p , for all x G M 2 ; 

(ii) is a regular poini0 of the spectrum of the linear operator H = — A + V acting on L 2 ; 

(iii) H acting on L 2 has exactly one negative eigenvalue Eq < with corresponding normalized eigenvector 
00- It is well known that ipa(x) can be chosen strictly positive and exponentially decaying as \x\ — > oo. 

Conditions (i)-(ii) guarantee the applicability of dispersive estimates of Murata |T7] and Schlag [55] to the 
Schrodinger group e~ lHt . These estimates are used for obtaining Theorems 14.11 and 14.21 see also [13j section 
4]. In particular (i) implies the local well posedness in H 1 of the initial value problem ljl.lll.2|l . see section [3] 
By the standard bifurcation argument in Banach spaces 19J for ()2.2() at E = Eo, condition (iii) guarantees 
existence of nontrivial solutions. Moreover, these solutions can be organized as a C 2 manifold (center manifold) 
for x G R", see [HI section 2] or [TT]. The proofs for the following results can be found in [TU section 2] or 

Proposition 2.1 There exist 8 > 0, the C 2 function 

h : {a G K x K : \a\ < 6} i-> L 2 n i/ 2 , a G K 

anrf t/ie C 1 function E : (—6,6) > K swc/i i/iai /or |J5 — £?o| < ^ KV'O) V'-b)! < <5 t/ie eigenvalue problem 
112. Sty has a unique solution up to multiplication with e %B , 8 G [0, 27r), which can be represented as a center 
manifold: 

ip E = aipo + h(a), E = E(\a\), (ipo, h(a)) = 0, /i(e ie a) = e ie /i(a), |a| < <5. (2.3) 

Moreover E(\a\) = C(|a| 1+Ql ), /i(a) = C(|a| 2+Ql ), and for a el, |a| < <5, /i(a) is a real valued function with 
0(a)=O(|ar) ondf (0) = 0. 

Since '00 (^) is exponentially decaying as |x| — > oo the proposition implies that 0b G L 2 ,. A regularity 
argument, see [25], gives a stronger result: 

Corollary 2.1 For any a € K, i/iere exists a finite constant C a such that: 

|| < X >" i\)e\\hi < Ccr\\lpB\\lP- 

We are going to decompose the solution of (|l.l[) - (|1.2p into a projection onto the center manifold and a 
correction. To insure that the correction disperses to infinity on long times we require that the correction is 
always in the invariant subspace of the linearized dynamics at the projection that complements the tangent 
space to the center manifold. A short description of the decomposition follows, for more details and the proofs 
see [T2]. 

Consider the linearization of (jl.ip at a function on the center manifold ipE = aipo + h(a), a = a% + ia2 G 
C, \a\<6: 

— = -iL^ E [w] - iEw (2.4) 

where 

Lip B [w] — (—A + V - E)w + Dg 4 , E [w] = (-A + V - E)w + lim + ew) - g($ E ) (2 g) 



1 see [23] Definition 7] or Af M = {0} in relation (3.1) in [17] 
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Remark 2.1 Note that for a G R we have we — &Wo + h(a) is real valued and 

Dg^ E [w] = 9 '{weWw + i^^he = \ [q'^e] + ^j^j w+ l - [g'{^ E ) 



— ; I W 

We 



hent 



\Dg^ E [w}\ < \w\max\ \g'(tp E )\, 



9 {we) 



Wi 



<C(\w E \ L+ai + \w E \ L+a2 )\M (2.6) 



where we used (jl.3p . For a = \a\e l9 G C we have, using the equivariant symmetry (|1.4|) . ip E — aipQ + h(a) = 
e z6 (\a\wo + h(\a\) = e ie ip E al , where ?/>^ al is real valued, and Dg^ E [w] = e %e Dg^t&,i[e~ %0 w], hence (|2.6p is valid 
for any ip E on the manifold of ground states. 

Properties of the linearized operator: 

1. L^ E is real linear and symmetric with respect to the real scalar product 3?(-, •), on L 2 (R 2 ), with domain 
iJ 2 (R 2 ). 

2. Zero is an e- value for —iL^ B and its generalized eigenspace includes i^ 5 -, } 

3. span K {|^-, an< ^ ^ a = { — *f^7' *l$rf} > wnere orthogonality is with respect to the real scalar 
product in L 2 (R 2 ), are invariant subspaces for —iL^ E and, by possible choosing S > smaller than the 
one in Proposition 1 2. 1[ we have: 

L 2 (R 2 ) = span R §^1 © H a , for all a G C, \a\ < 6. 

[ aa\ 002 J 

Note that Ho coincides with the subspace of L 2 associated to the continuous spectrum of the self-adjoint 
operator H = — A + V. 



4. the above decomposition can be extended to H 1 (R 2 ) 

dipE chpE 
dai ' da2 



H-\R 2 ) = span R \ ^ } © H a , for all a G C, \a\ < S, (2.7) 



where 



H a = U e H- 1 I 5R(-z^, <j>) = 0, and 3t(i^, <j>) = o) 
^ oa2 oa\ J 



Our goal is to decompose the solution of fll.l[) at each time into: 

u = We + V — a Wo + H a ) + n, n G 7Y a 

which insures that 77 is not in the non-decaying directions of the linearized equation (|2.4| at we ■ The fact that 
this can be done in an unique manner is a consequence of the following lemma: 

Lemma 2.1 There exists 5/2 > Si > such that any (f> G H^ 1 (R 2 ) satisfying W^Wh- 1 < ^1 can be uniquely 
decomposed: 

<P — We + V = a-Wo + h(a) + n 

where a = a\ + 102 G C, \a\ < <5, n G H a - Moreover the maps <f> 1— > a and <f> ^ r\ are C 1 and there exists the 
constant C independent on (p such that 

\a\<2U\\ H -i, \\tjWh-i <CU\\ h -i, 

while for 4> G L 2 (R 2 ) we have n G i 2 (R 2 ) and: 

\a\ < 2\\</>\\v, h|| L2 <C||0|| L2 . 
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Remark 2.2 The above lemma uses the implicit function theorem applied to 

F : R 2 x H-^M 2 ) ^ R 2 F(ax,a 2 ,<p) -- 
where ipE — (ax + ia 2 )ipo + h(a\ + ia 2 ) and 



da2 \ da 2 ' da\ 



-l 



da\ \ da\ da 2 

form the dual basis of j , | with respect to the decomposition (|2.7[) . ./Vote £/ia£ 

<9.f 

(ax, a 2 , 4>) = I R 2 - M(a x ,a 2 , 4>) 



d(ax,a 2 ) 

where the entries of the two by two matrix M are 



Mij = 9? (-57^, 



and, consequently, M(0,0, 0) is the zero matrix. Thus the implicit function theorem applies to F = 0, in a 
neighborhood of (ax, a 2 , (f>) — (0,0,0) and the number 5x in the above lemma is chosen such that: 



Vt(i d ^ E 



> — , whenever |(ai,a2)| < 2Sx, 



da\ ' da 2 

and the norm of the matrix M as a linear, bounded operator on R 2 satisfies: 

UUW = \\M(ax(<f>),a 2 (<f>),^\\ < \, whenever H^H^-x < Sx, (2.8) 

see )12\ section 2] for details. 

We need one more technical result relating the spaces 7i a and the space corresponding to the continuous 
spectrum of — A + V : 

Lemma 2.2 With 5\ given by the previous lemma we have that for any a € C, \a\ < 2Sx, the linear map 
Pc\n a '■ 7~La | — > 7~to is invertible, and its inverse R a ■ T~Co | — > 7~ia satisfies: 

||#aCIUl CT ^ C-4(\\l*_„, a e R and for all ( e Ho n L 2 _ a 
11-RaClli* < CpIICIIlp, 1 <P < oo and for all C G H n L p 

where the constants C- a , C p > are independent of a € C, \a\ < 25x- 
We are now ready to prove our main result. 

3 The Main Result 

Theorem 3.1 If hypothesis (fO]) . ([T3]) . (HI) hold and 

1 

2 <ai 

then there exists q' < ^^al ant ^ e o > such that for all initial conditions uq(x) satisfying 

max{||u || L ^, ||wo||i?i} < £o 
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the initial value problem hl.l}) - i ("Oj) is globally well-posed in H , and £/ie solution decomposes into a radiative 
part and a part that asymptotically converges to a ground state. 

More precisely, there exist a C 1 function a:Ri->C such that, for all t eM we have: 

u(t, x) = a(t)i{j n (x) + h(a(t)) +r/(t, x) 
v * ' 

where ipE(t) is on the central manifold (i.e it is a ground state) and r){t,x) G W a ( t ), see Provosition \2.1\ 
and Lemma \2.1\ Moreover, there exists the ground states ipE± oa o,nd the C 1 function 6 : R i— > R such that 
limiti^oc 9(t) = and: 

lim UE{t)~e- lt ^- e ^H E± J\ H 2r, L 2 = 0, <jeU (3.1) 



while r\ satisfies the following decay estimates. Fix po > max{ ai _ 1 ^ 2 > (4 + 2ct2) go -^4+202) J' 1 where qo 
-^-r > 4 + 2a 2 . Then for 2 < p < p ° qo „ we /iaue: 

1-2/p 



,0 l O g 1 - 2 /P0 (2+|t|) -r \ -1 -1 j 2 

Ceo (1+ u|)iI27p if ai > 1 or en < 1 and p < j^f 



(i+\t\y-i/p - Ui U1 ^ " ^ l ^ - i-«i+2/ P0 ' 

W*)IU» < <! (3.2) 

°l-2/p 

r log i- 2 ^ (g+]t|) - f , j 2 

° e ° (i+|t|)°i- 2 /po it ai < 1 and p > 1 _ Ql+2/po , 
/or some constant C — C(po). 

Remark 3.1 The estimates on rj show that the component of the solution that does not converge to a ground 
states disperses like the solution of the free Schrddinger equation except for a logarithmic correction in LP 
spaces for critical and supercritical regimes, ot\ > 1. In subcritical regimes, d\ < 1, the decay rate remains 
comparable to the free Schrddinger one in L p spaces for 2 < p < 2/(1 — a\), while it saturates to |t| Ql_1_0 in 
W , p > 2/(1 - ai ). 

Proof of Theorem 13.11 It is well known that under hypothesis (Hl)(i) the initial value problem (jl.ljl - 
(|1.2|) is locally well posed in the energy space H 1 and its L? norm is conserved, see for example Corollary 
4.3.3. at p. 92]. Global well posedness follows via energy estimates from HuoIIh 1 small, see Corollary 6.1.5 
at p. 165]. 

We choose £0 < b~\ given by Lemma [2.11 Then, for all times, ||u(£)||#-i < ||u(i)||x,2 < £0 < <5i and, via 
Lemma 12.11 we can decompose the solution into a solitary wave and a dispersive component: 



u(t) = a(t)ip + h(a(t)) + n(t) = if) E (t) + r)(t), where \a(t)\ = |ai(t) + ia 2 (t)\ < 2e Q < 28 x Vt G R. (3.3) 

Note that since h(a) is C 2 , see Proposition ^. 1[ and a is uniformly bounded in time we deduce that there 
exists the constant Ch > such that: 

m«|||^(*)||fl»J|^(<)||^,||^(*)||^}<C! ff eo J foraUiGM, 



which combined with Corollary 12.11 implies that for any a G R there exists a constant Ch,<j > such that: 

max (|| <x> a i) E (t)\\ H 2,\\ <x> a ^(t)|| H2 ,|| <x > a ^{t)\\m\ < C H .*e , for all t G R. (3.4) 
( oa\ 002 J 

Consequently, using the continuous imbedding H 2 (R 2 ) ^ L P (R 2 ), 2 < p < 00 and L 2 (R 2 ) ^ L : (R 2 ), a > 1 
we have that for all 1 < p < 00 and all a G R, there exists the constants C PiCT such that 



sup max 

ten 



H^f ll^^C*)!!^' ll*i(«(*))II^S» < Cp^eo, (3-5) 



see Remark |2. 21 for the definitions of ^j(a), j = 1,2. In addition, since 

ue C(R,i? 1 (R 2 ))nC 1 (R,i/" 1 (R 2 )), 
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and «ho respectively u i— > 77 are C 1 , see Lemma I^TTT we get that a(i) is C 1 and ry € C(R, i/ 1 ) n C 1 (R, H 1 ). 

The solution is now described by the C 1 function a:M^C and 77(f) € C(K, ff 1 ) n C^R, if -1 ). To obtain 
estimates for them it is useful to first remove their dominant phase. Consider the C 2 function: 

9(t) = f E(\a(s)\)ds (3.6) 
Jo 

and 

u(t) = e l9{t) u{t), (3.7) 

then u(t) satisfies the differential equation: 

idu(t) = -E(\a(t)\)u(t) + (-A + V)u{t) + g(u(t)), (3.8) 

see (jl.l j) and (|1.4p . Moreover, like u(f), can be decomposed: 

u(t) = o(t)V»o + fr(5(t)) +»)(*) (3.9) 
s . ' 

where 

~a(t) = e ie ^a(t), fj(t) = e^r,(t) e W 8(t) (3.10) 
By plugging in (|3.9[) into (|3 .8[) we get 

iJl +iDq p E \.^ = (-A + v-E(\a\)$ B + v)+g$E)+g$B + ®-g$E) 

or, equivalently, 

dry d^ E da x di> E da 2 . _ . , 7 , . 



Gs pan K {^,^f} 



where is defined by (|2.5|) : 



L ^ = (-A + V - £(|5|))»7 + ^- £ g^E + erT)| e=0 



and we used |a| = |a|, while g 2 is defined by: 



52(-0B,?y) = fl#s + »7) - - -fe.9{$E + £fi)\e=o (3.12) 



and we also used the fact that -0s is a solution of the eigenvalue problem (12. 2|) . Note that 172 is at least 
quadratic in the second variable, more precisely: 

Lemma 3.1 There exists a constant C > such that for all a, z G C we have: 

\g2(a,z)\<C(\ar + \ar + \zr + \z\ aa )\z\ 2 
Proof: From the definition (|3.12p of g 2 we have: 

g 2 (a, z) = g(a + z) - g(a) - Dg a [z] = / (Dg a+TZ - Dg a ) [z]dr = / / D 2 g a+STZ [Tz][z]dTds. 

Jo Jo Jo 

Now (ll.3[) and (| 1 .4[) imply that there exists a constant C\ > such that the bilinear form Dg on C x C 
satisfies: 

\\D 2 g b \\ <Ci(H Ql + |&r a ), V6eC. (3.13) 
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Hence 



\9 2 (a,z)\ < d ((2max(|a|, + (2max(|a|, \z\)) a *) -\z\\ 



which proves the lemma. □ 

We now project p. lip onto the invariant subspaces of — iL^ E , namely span M {^^-, ^j-}, and 7i a . More 
precisely, we evaluate both the left and right hand side of (|3.11j) which are functionals in H~ 1 (M. 2 ) at = 
^j(d(t)), j — 1,2, see Remark 12.21 and take the real parts. We obtain: 





d 


Si 








+ dt 


0,2 







where 

92j{^E,f]) = $l( i &j,-ig2(i>E, fj))> 3 = !, 2 - 
Note that from Lemma 13.11 and Holder inequality we have for all t G M. : 



(3.14) 



\92 3 {^E{t),m)\<c [ \y 3 (t,x)\(\^ E (t, X )r + \j E (t 7 x)r* + \f!(t,x)r + \f!(t 7 x)rA m, X )\ 2 d X (3.15) 



c 



ll*j(*)IUt, II^WII^ + II^WIIl- 11^)111- + ll*i(*)IU-i|l^)ll2t 2 ai + 11^)11^11^)1 



|2+Q 2 
ILP2 



where r^ 1 + {p 2 /2)- 1 = 1, r^ 1 + (p 2 /(2 + c^))" 1 = 1, 3 = 1, 2. 

To calculate 3?(\l>j, ||), j = 1,2 we use the fact that 7j(£) S 7Y a , for all fel, i.e. 

Differentiating the latter with respect to t we get: 



which replaced into above leads to: 

d 
dt 



S 2 



(I R2 - Ms)" 1 



92\{ilE,ri) 

g22(-4>E,v) 



(3.16) 



where the two by two matrix M„ is defined in Remark 12.21 In particular 

= -M fi (I R2 - M^)- 1 



922{i>E,f)) 



which we use to obtain the component in H a = {*i(a), \E , 2(S)}" L of (|3. 1 If) : 

-qT = -iL$ B fj-ig2(i>E,fj) - (I-Mu)~ g3(i>E,fj), 
where 173 is the projection of — ig<z onto span R {-^^-, relative to the decomposition (|2.7I 

93(-4>E,v) = 92l(lpE,v)^^- +522(^B,??)^J, 



(3.17) 



see (|3.14p for the definitions of g^j, j — 1,2, and I — M^ is the linear operator on the two dimensional real 
vector space span R {^^, ^J-} whose matrix representation relative to the basis {^^-, ^7} is Irz — M„. It 
is easier to switch back to the variable rj(t) = e~ ie ^fj(t) e TL a : 



— = -i(-A + V)rj - iDg^ B rj ~ ig 2 {'>pE,i) - (I - M u )~ g 3 (t/) E , 77), 



(3.18) 
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where we used the equivariant symmetry (|1.4p and its obvious consequences for the symmetries of Dg, g 2 -, g$ 
and M. Since by Lemma l2~2l it is sufficient to get estimates for ((t) — P c rj(t), we now project (|3.18|) onto the 
continuous spectrum of — A + V : 

dC 

-± = -i(-A + V)( - iP c D gi , E RaC - iPc92(i>E, RaC) - PS - M u )- x gz{i> E , #aC), (3-19) 

where R a ■ "Ho > Ti-a is the inverse of P c restricted to Tt a , see Lemma \2 .21 
Consider the initial value problem for the linear part of (|3.19[) : 

o z 

— = -i(-A + V)z-iP c Dg^ B{t) R a(t) z(t) (3.20) 
z(s) = v e Ho 

and write its solution in terms of a family of operators: 

Q(t,s) :Ho^Hq, Sl(t,s)v = z(t), t, s € K. (3.21) 

In Section U we show that such a family of operators exists, is uniformly bounded in t, s with respect to the 
L 2 norm and it has very similar properties with the unitary group of operators e~ l (~ A+v ^ t ~~ s * > P c generated 
by the Schrodinger operator — i{— A + V)P C . In particular fl(t, s) satisfies certain dispersive decay estimates 
in weighted L 2 spaces and L p , p > 2 spaces, see Theorem 14.11 and Theorem 14.21 For all these results to hold 
we only need to choose £o small enough such that £oCh,4o-/3 < ei, where er > 1 and £\ > are fixed in Section 
4 and the constant Chau/3 is the one from (|3.4p . 

Using Duhamel formula, the solution C € C(E, H 1 n W ) H C* 1 ^, ^(R 2 ) H H ) of (j3~Tgj) also satisfies: 

C(t) = n(*,0)C(0)-i f n(t,s)P c g 2 (ip E (s),R a ( s ){(s))ds 
Jo 

Q(t, s)P c (I - M„( s ) ) _1 p3 (V'b (s) , E a(s) C(s))ds. (3.22) 

Note that the right hand side of (|3 . 22[) contains only terms that are quadratic and higher order in (, see 
Lemma |3. II and (|3.15|) . As in [131 112] this is essential in controlling low power nonlinearities and it is the main 
difference between our approach and the existing literature on asymptotic stability of coherent structures for 
dispersive nonlinear equations, see [T31 p. 449] for a more detailed discussion. 

To obtain estimates for £ we apply a contraction mapping argument to the fixed point problem (|3.22p in 
the following Banach space. Fix po > 2 such that 

p >max{ 2 (4 + 2a 2 ) gC ~ 2 - } , (3.23) 

Lai -1/2 <7o-(4 + 2a 2 )J 

and let 

P2 = : 3 - 24 

Po + qo - 2 

and 

Pl = - 2 -^, ifai<l, (3.25) 

1 - ai + 2/po 

then 

Case I if «i > 1, or 1/2 < a% < 1 and pi > let: 

y= ^£C(l,L 2 nl ft ) : Bup||«(t)||ia <oo, sup (1 + ^ *\ ||«(t)|Up2 < c 

[log(2+|i|)] 1 -™ 
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Case II if 1/2 < ai < 1 and p\ < p 2 , let: 



Y = IveC (R, L 2 n L Pl n L P2 ) : suplK*)!!^ < oo, 

tGR 



(i + KI)^" II /.Ml 

sup ^H^'WII^ 1 < 00 > SU P' 

tGR pi tGR 



(1 + 1*1 



|w(t)||i)»a < OO 



[log(2 + |t|)P 



pog(2 + |t|)3 l " 



endowed with the norm 



(1 + 

\v\\ Y = max •( sup ||w(t)|| L 2, sup —^\\v{i)\\LP2 

tGR tGR P2 



[log(2 + |t| 



in Case I, while in Case II 



\v\\y = max < 



sup\\v(t)\\ L2 , sup {1 + m> . 11^)11^, sup 1 +l U - _^ H*)ll^ 



/gp 



I^+ltl)] 1 - 
Consider now the nonlinear operator in (|3.22D : 



[log(2 + |t|)] l "5o 



N(v)(t) = -i Q(t,s)P c g 2 (Ms),R a {s)v(s))ds- Q(t,s)P c (l - M^y 1 g 3 (ip E (s), R a ( s ) v ( s ))ds. (3.26) 
Jo Jo 

We have: 

Lemma 3.2 N : Y — > Y is well defined, and locally Lipschitz, i.e. there exists C > 0, such that 

\\Nv x - Nv 2 \\ Y < C(\\ Vl \\ Y + \\v 2 \\ Y + |M|^ +Ql + + IM| y +Q2 + ||t' 2 ||^ +a2 )lki - v 2 \\ Y . 

Assuming that the lemma has been proven then we can apply the contraction principle for (|3.22p in a 
closed ball in the Banach space Y in the following way. Let 

« = n(t,o)c(o) 

then by Theorem 14.21 

||«|| y < max{C2,Cp o , Pl ,Cp o ,p 2 }||C(0)|| L2nL ^ 

where we used the interpolation ||C(0)l!i r < l!C(0)|| i2 L i' i Qo — r — 2 with r = q' and r = p 1 defined in 
theorem 14.21 for p = pj, j = 1,2. Recall that 

C(0) = P cV (0) = P c u - %(0)) =u - {i/>o,U )i> - h(a(0)) 

where uq = u(0) is the initial data, see also (|3.3p . Hence 

||C(0)|| £an£ ^, < H«o|| £an£ «&+||«olU>ll^ll £ ^ +£>i||«o|U= < De 

where D\, D > are constants independent on uo and the estimate on h(a(0)) follows from Proposition 12. II 
and |a(0)| < 2||m ||l 2 see Lemma |2~T1 

Therefore we can choose So small enough such that R = 2\\v\\y satisfies 

Lip d = f 2C(R + R 1+ai + R 1+a2 ) < 1. 
In this case the integral operator given by the right hand side of (13. 22[) : 

K(0=v + N(0 
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leaves i?(0, R) = C £ Y : \\(\\y < R invariant and it is a contraction on it with Lipschitz constant Lip defined 
above. Consequently the equation (|3.22p has a unique solution in B(0,R) and because £ C(R,H 1 ) <^-> 
C(R, L 2 , L Pl , L P2 ) already verified the equation we deduce that ((t) is in B(0, R), in particular it satisfies the 
estimates (|3 . 2[) . 

Then r](t) = R a (t)((t) satisfies (|3.2[) because of Lemma [2.21 Moreover, the system of ODE's (|3.16[) has 
integrable in time right hand side because the matrix has norm bounded by 2, see (|2.8p . while g 2 j satisfy (|3.15[) 
where fj(t) differs from rj{t) by only a phase and the L p , 1 < p < oo norms of ^/j(t), ipE(t) ar e uniformly 
bounded in time, see (|3.5p . Consequently cii(t) and a 2 (t) converge as t — > ±oo, and there exist the constants 
C, e > such that: 

lim o(t) = lim aAt) + ia 2 (t) = a ±00 , |a(±i) - a ±oc | < C(l + i)~ £ , for alii > 0. 

t — >±oo t — >±oo 

We can now define 

fe ±oc = a±oo-0o + h(a ±00 ), (3.27) 

and we have 

lim \\<4> E {t) - ^ E±oa \\h>m.' = 0, for a £ R (3.28) 

t— >±oo 17 

where we used (|3.9[) and the continuity of h(a), see Proposition l2.il In addition, since E : \—28\, Si] i— > (—5, 5) 
is a C 1 function, see Proposition ^. li the following limits exist together with the constant C\ > such that: 

lim E(\a(t)\) = E ±oa , \E(\a(±t)\) - E ±oa \ < d(l + for alH > 0. 

£ — >±oo 



If we define 



then 



± J* E(\a(s)\) -E +00 ds if t > 
0(t) = <( if t = (3.29) 

A /*£(|a(fl)|) if t < 



lim 6>(t) = 



and 



w =/W)ik, ={;<£ : :*™ (3.30 

where we used \a(t)\ = \a(t)\, see (|3.10j) . 

In conclusion, since ip E (t) = e i6 ^^ E (t), see (031), (JSH]) and ([3~TU]) , we get from (pH5|) and ([330)) the 
convergence (|3.1[) . 

It remains to prove Lemma 13.21 

Proof of Lemma [372} It suffices to prove the estimate: 

\\Nvi - Nv 2 \\ Y < C(|M|y + |M|y + |M| r +ai + lk 2 ||^ +Ql + IMI y +a2 + lkl| y +Q2 )||^i - v 2 \\ Y , (3.31) 

because plugging in v 2 = and using AT(0) = 0, see (|3.26l) . will then imply N(vi) € K whenever i>i e K 
Note that via interpolation in L p spaces we have for all v £ Y and any 2 < p < p 2 : 

l-2/p 

I, log !- 2 /po (2+|t|) -r \ i ^ i i ^ 2 

M|y (i+| f |)i-a/P if ai > 1 or ai < 1 and p < 
K*)llz» < <! (3-32) 

oil-2/po 

|| log 1-2/PO (2+|t|) -r . 1 , 2 

Nlr (i + | t |)°i-2/po ifa 1 <landp> 1 _ Q1+2/po - 

Now, from (|3.12p , we have for any V\, v 2 £ Y : 

g 2 (ipE,R a vi) - g 2 (ipE,R a v 2 ) = g(ipE + R a v\) - g(ip E + R a v 2 ) - Dg^ E [R a {vi - v 2 )] 

1 

(Dg^E+Rairvi+ii-T^) - Dg^z) [R a (vi - v 2 )}dr 
1 ,1 

/ D 2 9-<PE+sB.a(Tv 1 +(i^T)v 2 )[Ra{Tv l + (l - t )v 2 )] [R a (vi - v 2 )]drds 
Jo 

Ai(ipE,vi,v 2 ) + A 2 (tp E ,vi,v 2 ) + A 3 (ipE, vi,v 2 ), (3.33) 
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where we consider Xj{t, x ), 3 = 1, 2 to be the characteristic function of the set S\ — {(t,x) G 

\ipE(t,x)\ > max(|i? a(t) wi(t,a;)|, \R a(t) v 2 (t,x)\)}, respectively 

S 2 = {(t,x) eRxR 2 : iiia3t(|J? o(t )t;i(t,a:)|,|i2 o(t) tj2(t,ar)|) < 1} and 



x 



>-2 



1 pl 

2, 



A 1 (ip E ,v 1 ,v 2 ) = Xi D 9i,E+sRa.(Tv 1 +(i-T)v 2 )[ R a(Tv 1 + (1 - T)v 2 )][R a (vi - v 2 )]drds, 

Jo Jo 

,vi,v 2 ) ~ (1 — Xi)X2-D g^ E + S R a (Tv 1 +(i-T)v 2 )[RaTVi + (1 — T)v 2 )][R a (vi — v 2 )]drds, 



'0 Jo 

A 3 (tp E ,ui,u 2 ) = (l-Xi)( 1 -X2)D 2 gip E+S Fi, a ( TVl+ (i- T ) V2 ){Ra{TVi + (l-T)v 2 )][R a (vi-v 2 )}dTds. 

Jo Jo 

Note that there exists a constant C > such that for any ip E , v\, v 2 G Y, any tgl and almost all x € M 2 
we have the pointwise estimates: 

\Ai(ip E (t,x),vi(t,x),V2(t,x))\ < C(2 a ^ E (t,x)\ ai +2 a2 \^ E (t,x)D(\R ait) v 1 (t,x)\ + \R a{t) v 2 (t,x)\) 

x\R a {t)(vi{t,x) - v 2 (t,x))\ 

\A 2 (ip E (t,x),Mt,x),v 2 (t,x))\ < 2^C{\R a{t) v 1 (t 1 x)\ 1+ ^ + \R a{t) v 2 (t,x)\ 1+ai ) \R a{t) ( Vl (t,x) - v 2 (t,x))\ 
\A^ E {t,x), Vl {t,x),v 2 {t,x))\ < 2 a -C(\R a(t) v 1 (t,x)\ 1+a2 + \R a (t)V 2 (t,x)\ 1+a2 ) \R a{t) ( Vl (t,x) - v 2 (t,x))\ 

where we used (|3 . 1 3|) . Consequently, for any ueM there exists a constant C a > such that for any tgR: 

11^1(^(0, «i(*),«2(*))IU» < C||2 Ql |^W| ai +2 a2 |VBWI a2 IU;(ll«a( t )«i(*)IU~ +\\R a{ tMt)\\L*2) 

x\\R a{t) (vi(t)-v 2 (t))\\ LP2 

^ Ctr (°i g +|t])t QHI y + - Mr (3-34) 

(3.35) 

where 7 + ^ = 5; an< 4, for j = 1,2 defined in Remark 12.21 

\&(*Mt)),-Mi(Mt)Mt)Mt)))\ < \\^Mm\LiJAi(Mt),vi(t),v 2 (t))\\Li 

CAog ai (2 + \t\) ru 
- 62 (1 + u|)b! (\\vi\\y + \\v 2 \\y)\\vi - v 2 \\ Y (3.36) 

where 

f 2-i in Case I, f 2±^§^ in Case I, 

fc-U = i*C-.n: inCase I, 

see the definition of the Banach space Y, and we used Holder inequality together with (|3.5[) and Lemma 12.21 
Similarly, for any 1 < r' < 2 we have (2 + ai)r' < (2 + o^V < P2, hence the above pointwise estimates 
and (13.321) imply that there exists a constant C r > > such that for any t £ R : 

\\A 2 ^ E (i),vi{t),v 2 {t))\\ Lrl < 2 a >C\\ |i? o(t)Ul (i)| 1+ai + |i? o(t) t; 2 (i)| 1+ai || (3 + aiV ||iZ, (f) f^(t) -«2(t))ll r , (a+ a l)r / 
where 

2 „ _ ai+2/r 



& 2 (r') = ai + a 2 (r') = 2i±^I, if a x > 1 or a x < 1 and (2 + ai )r' < Pl , 

b 2 (r>) = (2 + a 1 )(a 1 -±), a 2 (r') = (2 + ai )^^-, if «i < 1 and (2 + ai)r' > Pl , 



with 1/?' + 1/r 1 = 1, and 

\\M$E®Mt)Mt))\\Lr> < 2 a2 C|| |i?a W ^lW| 1+Q2 +l^a( 4 )«2W| 1+Q2 || \\R a(t) (Vl (t) - V 2 (t)) || r , (a+tta )r , 



L i+»2 



^ Cr (H|li)£ - (H^llr" 2 + WM^vi -Mr, (3.40) 
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where 

b 3 {r') =a 2 + §, 



b 3 (r') = (2 + a 2 )( ai -±), a3 (r') = (2 + a 2 )^^, if ax < 1 and (2 + a 2 )r' > Pl 



if Q!i > 1 or oli < 1 and (2 + a 2 )r' < pi, 

at— 2/po • r , , ,„ . w : (3-41) 



Moreover, using Cauchy-Schwartz inequality and (|3.5[) we have: 

C 2 log a ^ 2 >(2+ |t|) /M ^ „!+,,,, 



< c 2 . 



(l + \t\) b ^) 



-(\\vi\\? ax + \\v 2 \\? ai ) \\vt -v 2 H0.42) 



and 



|jR(^(a(t)),-a 3 (^(t),t;i(A)^2(t)))l < C 2 ,o 2 „ g , u ,l^ (IHll^ 2 + IHl^ + ° 2 )lki - v 2 \\ Y . (3.43) 



(l+|i|) b 3(2) 



Now, from (|3~T7|) and (pTH)) we have 

g 3 (4>E, R a vi) - ga(il>E, R a v 2 ) 



dib dip 

= ^{^!i{a),-i(g 2 {^E,RaVi) - g 2 (ipE, R a v 2 )))^ h U(^ 2 (a), -i(g 2 (ipE, R a vi) - g 2 {tps, R a v 2 )))-^ — 

oai oa 2 

= R<*i(o), -i(Ai + A 2 + A 3 )(^ E , Vl ,v 2 )) d -p- + »<*a(o), + A 2 + A 3 ){^ E ,v u v 2 )))^. 

oai oa 2 



Consequently, for 



A4(ipE,vi,v 2 ) = (I - M u ) 1 {g 3 {^E,R a vi) - g 3 (ipE,R a v 2 )) 



(3.44) 



we have that for any a £ K there exists a constant Co- > such that: 



\A4(<ip E {t),Vi(t),v 2 {t))\\ L 3 < max 



dai 



Ll: 



da 2 



< 



x^/\^ 1 (a(t)),-z(A 1 +A 2 + A 3 )(tW + |»(* 2 (a(*)),-i(^i + ^ + A 3 )(t))\ 2 
C CT log a4 (2 + \t\) nK _ , „__ , |U ||1+Q1 .. „i +a!1 „ |,i+ Q2 , n , _Mi+o 3 



(1 



,„ (||«i||r + IM|y + \\vi\\ Y +ai + \\v 2 \\^ + |H|^ Q2 + \\v 2 \\^)\\ Vl - v 2 \\ Y (3A5) 

(3.46) 



where 

b 4 = min{6i, 6 2 (2), 6 3 (2)}, a 4 = max{ai, a 2 (2), a 3 (3)}, 

and we used (j2l3j) . ([336]) . ([3^2]) . and ([3~43f . 

We are now ready to prove the Lipschitz estimate for the nonlinear operator N, (|3.31|) . From its definition 
(|3.26p and (|3.33[) . (|3.44p we have for any v\, v 2 € Y, any 2 < p < p 2 , and a fixed a > 1 : 



Q{t, s)P c (-«Ai - iA 2 - iA 3 - A 4 )(ip E (s), vx(s), v 2 (s))ds 



L'< 



< 



||^(M)IU=^i* {\\A 1 ^j e (s),v 1 ( S ),v 2 ( S ))\\ L 2+\\A 4 (^e(s),v 1 ( S ),v 2 (s) 



\Ll 



ds 



+ / \\tt(t,s)\\ Lq > nLP >„ LP (\\A 2 (iPE(s),v 1 (s),v 2 (s))\\ Lq ' nLP ' + \\A 3 (^ E {s),v 1 {s),v 2 (s))\\ Lq ' nLp/ )ds. 



where 

1/p' + l/p = 1, q' = p'(po - 2)/(po - p'), V<7 + 1/9' = 1- 
From Theorem 14.11 and estimates (|3.34[) . (|3.45|) we get: 

1*1 

\\n{t,s)\\ L 2^ LP (\\A 1 (l/} E (8),V 1 (8),V2(8))\\ L l + \\A 4 (ijE(s),v 1 {s),v 2 {s))\\ L 2)ds 
< (\\vi\\y + \\v 2 \\y + \\vi\\ Y +ai + \\v 2 \\ Y +ai + |h|| r +a2 + IM| y +a2 )|K -v 2 \\ Y 

ds 



(3.47) 



|t- S |l-2/P 



aio g ai (2 + |s|) aio g a4 (2 + |s| 



(i + N)^ 



(1 + H)<* 
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while from Theorem l4.2l and estimates (|3.38[) . (|3.40p we get: 

1*1 i i 

C po , p log^(2 + \t-s\) f^log Q2 ^(2+| S |) CVlog a2 ^(2+| S |)' 



| t _ s |l-2/p - (! + | s |)6 2 (g') ' (1 + \ s \)Up') 



ds 



and 



l-2/p 



C Po , p \og^(2 + \t-s\) f c q . log a3 ^ ) (2 + |s|) Cp> log a3(p,) (2+| S |) ' 

|f_ a |l-2/p maX 1 (1 + \ s \)b 3 W) ' (1 + | S |)Mp') 



ds. 



In Case I, i.e. ax > 1, or 1/2 < ai < 1 and pi > p 2 , since a 2 > «i and p 2 > 4 + 2a 2 > 4, we have from 
(13371) . f3~39l) . (13411 and f3~46)l for r' G {q',p',2} and 1/r+l/r' = 1 : 

h = 2 > 1, 6 2 (r') = «! + -> 1, 6 3 (r') = a 2 + - > 1, b 4 = min{6i, 6 2 (2), b 3 (2)} > 1. 

P2 r q 

We now use the following known convolution estimate: 

log n (2 + |f- a |) 1(^(2+1^ „ log a (2+|t|) 

(i + H) d " (a ' ' c ' j (i + |<|) h ' 

to bound the integral terms above and obtain for all 2 < p < p 2 : 



for d > 1, 6 < 1, (3.48) 



l-2/p 

logi-Vpo (2 4- tD 
||iVK) (t) N(v 2 )(t)\\ LP < C p \ l + ltl ]J 2 / p l) 

x(\\vi\\y + \\v 2 \\y + \\vi\\Y +ai + \\v 2 \\ Y +ai + \\vi\\ Y +a2 + |M| y +Q2 )||«i - %||y (3.49) 

(3.50) 

which, upon moving the time dependent terms to the left hand side and taking supremum over f £ 1 when 
p € {2,p 2 }, leads to (|3.31|) for C = max{C 2 , C P2 }. 

In Case II, i.e. 1/2 < ot\ < 1 and pi < p 2 , we have from (|3.37|) fei = 2(«i — — ) > 1 because p > 
2/(ai - 1/2), see ((3723|) . From ([339| . under the restriction 2 < p < p 1: with p', q' , q defined by ([347)) , we 
have either: 

6 2 (p') > 6 2 (g') = ai+2/ g > 1, 

or 

& 2 (p') = b 2 {q') = (2 + a x )(ai - 2/po) > (2 + e*i)/2 > 1. 
Since a 2 > cti implies &3(-) > & 2 (-) we deduce that, under the restriction 2 < p <pi, we also have 

fosb') > W) > W) > 1, 

and 

&4 = min{6i,6a(2),6 3 (2)}>l. 

We can again apply (13.48)1 to the above integral terms and get for 2 < p < p\ the estimate (13.491) . For p > pi 
one can show that (2 + a\)q' < p± hence b 2 (q r ) = ax + 2/g, and, in the particular case of p = p 2 , we get 

b 2 (q' 2 ) = ax + 2/q 2 < 1, 

where q 2 , </ 2 are given by (|3.47|) . We now have from convolution estimates: 

ll 1 — 2/p2 / ' \ 1 — 2/p2 | / m ct\—1 1 pq 

/•'*' log T ^(2+K- S |)log a2( ^ ) (2 + | S |) logT^+^fe)( 2 + | t |) log^^r( 2+ |t|) 

i |i-s|l-2/ P 2 (l + |s|)W S " lP2j (1 + |t|)ai+2/92-2/p a - ^ ( X + | f | )ai -2/ M ' 
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where we used (|3.47p and p 2 < Pq to obtain: 

2 2 _ 2 / l-2/p 2 \ k 2_ 

P2 Q2 PQ \1 - 2/poJ ~ p ' 

Since b 2 {p' 2 ) > 62(92) an d ^3(^2) ^ ^3(92) ^ ^2(92) we deduce 

ai-2/pp 
- log 1-2/fO (2 + \t\) 

\\N( Vl )(t)-N(v 2 )(t)\\ LP2 <C P2 ^^J/J* 

X(\\VI\\ Y + \\V2\\Y + WviW^ 1 + |k 2 ||y +Ql + IMIy^ 2 + ||«2 || y + ° 2 ) - «a||y 

(3.51) 

which, combined with (|3.49|) for p € {2,pi}, after moving the time dependent terms on the left hand side and 
taking supremum over t £ R, gives (|3.31|) in the Case II with C = max{C2, C Pl , C P2 }. 
This finishes the proof of Lemma 13.21 and of Theorem 13.11 □ 



4 Linear Estimates 

Consider the linear Schrodinger equation with a potential in two space dimensions: 

(i^ = (-A + V(x))u 
\u(0) = u . 

It is known that if V satisfies hypothesis (Hl)(i) and (ii) then the radiative part of the solution, i.e. its 
projection onto the continuous spectrum of H = —A + V, satisfies the estimates: 

\\e- lHt P c u4 L 2 <C M K \\uoh', S€l, (4.1) 

" °" L - " M (l + |f|)log 2 (2+|i|)" 1 ' 

for any a > 1 and some constant Cm > depending only on a see [TTJ Theorem 7.6 and Example 7.8], and 

\\e- lHt P c u \\ L v<^^\\u \\ LP , (4.2) 

for some constant C p > depending only on p > 2 and p' given by p'^ 1 + p^ 1 = 1. The case p = 00 in (|4.2|) 
is proven in [53]. The conservation of the L 2 norm, see [SJ Corollary 4.3.3], gives the p = 2 case: 

\\e- lHt P c uo\\ L i = \\u \\ L 2. 

The general result (|4.2p follows from Riesz-Thorin interpolation. 

We would like to extend these estimates to the linearized dynamics around the center manifold. In other 
words we consider the linear equation (|3.20p . with initial data at time s : 

Oz 

■g^ = - , i(-A + V)z-iP c Dg,p E{t) R a{t) z(t) 
z(s) = v E Ho 

Note that this is a nonautonomous problem as the bound state ipE around which we linearize may change 
with time. 

By Duhamel's principle we have: 

z{t) = e- lH ^P c v^i f e- lH(t -^P c D 9i , E(T) R a(T) z{T)dT (4.3) 

J s 

As in (|3.2ip we denote 

Q(t,s)v d = z(t). (4.4) 
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Relying on the fact that i/je (t) is small and localized uniformly in t G R, we have shown in [T31 Section 4] for 
the particular case of cubic nonlinearity, g(s) = s 3 , set, that estimates of type (|4.1|l - (|4.2|) can be extended 
to the operator Q(t, s). Due to (|2.6p which implies for a > and 1 < p' < 2 : 

H^ii.zlUj < c(||V B ||it Q1 +|Mit aa ) C_ CT ||z|| L2 (4.5) 

\ 2er/(l + ori) 2cr/(l + a 2 )/ 

il+ai | II < nl+a 2 » ~ " " ' ' 



\\Dg^ E R a z\\ LP , < C ll^li;!^,, +ll^li;| 1 ? a2 ), ) C^\\z\\ L 2 - = - + ^ (4.6) 
\\Dg^ E R a z\\ LP , < C(U E \\l+% l)q + \\^ E \\l+% 2)q ) C r \\z\\ L r, ^ = ~ q + l (4.7) 

see also Lemma [221 we can use, with obvious modifications, the arguments in [T3l Section 4] to show that: 

Theorem 4.1 Fix a > 1. There exists e\ > such that if || < x > 4<T / 3 ?/; s (£)|| ff 2 < ei for all t £ R, i/ien 
t/iere ea;is< constants C, C p > wii/i £/ie property that for any t, s G K i/ie following hold: 

W,S)\\ L 2^ L 2_ < 



and, for: 



(l + |i- S |)log"(2 + |i- S |) 
C 

Iflft, s)|| r p ' l _ > r2 < — 3-, for any 2 < p < oo where p' _ 1 + p~ 1 = 1 , 

{t-s^p 

\m,s)\\ Ll ^ LP < ^-t, for any p> 2 (4.8) 

|t-s| 1_ 5 

T(t,s) = n(t,s)-e- iH ( t -^P c , (4.9) 



Lemma 4.1 Assume that \\ < x > 4cr / 3 i/j E (t)\\jj2 < £i, t G K, where e\ is the one used in Theorem \4-l\ Then 
for each 1 < q' < 2 and 2 < p < oo i/iere exist the constants C q > , Cp,g' > such that for all t, s G K we ftawe: 



l r (*i s )|lLi n M'^i,i„ < 



i+i*-*r 

C Pig > log(2 + |t- 8 |) 
(l + lt-sD'-i 



Note that according to the proofs in [13l Section 4] C q i — + oo as g' ^ 1 and C Pi9 ' — > oo as o' — > 1 or p — * oo. 
These could be prevented and an estimate of the type 

IITtM)!!^- ^ 1 " 6 ' 2 ,:": 3 " (4.10) 

1 + \t — S\ 

can be obtained by avoiding the singularity of \\e~ lHt \\i J i^Loc 1 ~ t^ 1 at t = via a generalized Fourier 
multiplier technique developed in [12l Appendix and Section 4]. We choose not to use it here because it 
requires stronger restrictions on the potential V(x) like its Fourier transform should be in L 1 while its gradient 
should be in L p , for some p > 2, and convergent to zero as |x| — > oo. 

We now present an improved L? estimate for the family of operators T{t, s) : 

Lemma 4.2 Assume that \\ < x > 4<J / 3 iPe^Wh 2 < £ i> t G K, where £\ is the one used in Theorem \4-l\ Then 
there exists the constants Ci > such that for all t, s G M we have: 

\\T(t, S )\\ L 2„ L 2 <C 2 

Proof: We are going to use a Kato type smoothing estimate: 

|| < x e-^/CaOIUjp^) < C x ||/|| za , (4.11) 
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see for example [3T]. We claim that the previous estimate still holds if we replace e tH ^ S ^P C by Q(t,s), 
namely, there exists a constant Ck > such that for any s€l: 

|| < x >-° fi(-, s)f\\ Lm<Li) < C K \\f\\ L2 . (4.12) 

Indeed, from (|4.4|) and (|4.3p . we have 



< x >~ u fi(t, s)v =< x >-° e- H[t - s) P c v + J <x >-° e- tH{t - T) P c Dg 4 , E{T) [R a (T)n(T, s)v}dr 
and using (|4.5|) : 

< || e -«r(t— )„| Ua _,_ cr slip (||^(rr)||i+7 , + ll^)lli~7 fl+ ) 

TGl 2cr /( 1 + n l) 2(t/(1+q 2 )/ 

||fi(r, S )«( S )|U. 



(l + |i-r|)log 2 (2+|t-r|) dT - 
By Young inequality: ||/*3lU 2 (R) < II/IIlHR) IMIl 2 (R) and (|4.11j) we get 

s)w|| L 2 (RjL 2_ ct) < C K ||u|| L 2 + Cei\\il(; s)w|| L 2 (RjL 2 j 

which implies (|4.12p . 

Finally we turn to the estimate in L 2 for T(t, s) : 



\\T(t,s)v" 2 



= (f e- lH ^P c D gi , E [R a n(T, S )v]dT, f e^^-^PcDg^iRMiT^s^dT') 

J s J s 

= j j dTdT'(Dg^ E [R a n(T,s)v],e- iH( - T - T ')p c Dg^ B [R a n(T> »«]) 
< C S u P (||^(r)||i|f; +Il^(r)||it7 J 

r ^j^ \ 2a-/(l + a 1 ) 2a-/ (l + a 2 ) / 

x [ [ drdr'\\Q(T, S )v\\ L2 \\ e - lH ^ P C \\ L ,_ ||0(r', a )«|| L2 . 



SL 2 (R) GL 1 (R) ei 2 (R) 

Using (|4.ip combined with Young then Holder inequalities the integral above is bounded by 

Cm\\^(;s)v\\ 2 l2{m>l2 _j < C M C 2 K \\v\\ L 2. 

where, for the last inequality we employed (|4. 12[) . Consequently, there exist a constant C 2 such that for any 
t, s£l: 

\\T(t, S )v\\ L 2 < C 2 \\v\\ L 2. 

This finishes the proof of the Lemma. □ 

Fix now 2 < p < oo and let p' — Po/(Po — !)• By applying Riesz-Thorin interpolations to the operators 
T(t, s) satisfying for all t, s € R : 

\\T(t, S )\\ L 2^ L 2 < C 2 

\\T(ts)\\ , < Cpo log(2 + |t- S |) 

M(M)II LW _ LP0 _ (1 + |t _ jf|)1 _jL 

we obtain that for any 2 < p < p there exists a constant C P0jP between C 2 and C po such that: 

1-2/p 

wrpu mi ^ Cpo.p log 1 - 2 /^ (2 + |t-s|) , p , Po - 2 

l|T(M)ll^ n ^^ < (1 + |t _ g|)1 -| ■ where p = — p fl = P^7' 

Finally, using (|4.9|) and the estimates for the Schrodinger group (|4.2p we get: 
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Theorem 4.2 Fix 2 < po < oo and assume that || < x > 4cr / 3 ip E (t)\\ H 2 < ei, i 6 R where e\ is the constant 
obtained in Theorem \4-l\ Then there exists the constants C2, C POlP > such that for all 2 < p < po and 
t, s G K the following estimates hold: 

\\n(t, s )\\ L ^ L2 < c 2 - 

no/, mi . C po , p log(2 + |t-s|)^fe , p , po-2 

\Mt,s)\\ Lg > nLP >^ LP < ■ —2 , where p = — — , g = p 

t — s * " PO P 

Note that the estimates for the family of operators f2(t, s) given by the above theorem are similar to the 
standard LP 1— ► i p estimates for Schrodinger operators (|4.2p except for the logarithmic correction and a 
smaller domain of definition L q P\ L p C £ p where q' < p' when p' < 2. If we would have proven (|4.10p then 
we could use po = 00, hence g' = p' in the above theorem and obtain: 

won mi <r CpiogCa + lt-al) 1 - 2 ^ . , P 

ll fi (^s)llLP'^Lp < — : "TTa where p = -. 

It — s\ p P 1 
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